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Abstract 

A numerical study of static, spherically symmetric monopole solutions of a sponta- 
neously broken SU(2) gauge theory coupled to a dilaton field is presented. Regular 
solutions seem to exist only up a maximal value of the dilaton coupling. In addition 
to the generalization of the 't Hooft-Polyakov monopole a discrete family of regular 
solutions is found, corresponding to radial excitations absent in the theory without 
dilaton. 
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The aim of this paper is to present a detailed numerical study of classical solutions of an SU(2) 
Yang-Mills-Higgs (YMH) theory coupled to a (massless) dilaton (DYMH) , when the Higgs field 
is in the adjoint representation. Dilaton fields appear naturally in low energy effective field 
theories derived from superstring models |2|, ^ . As previous studies have already shown |Q, ^ 
the inclusion of a dilaton in a pure Yang-Mills (YM) theory has drastic consequences already 
at the classical level. In particular the dilaton Yang-Mills (DYM) theory possesses finite energy 
'particle-like' solutions which are absent in the pure YM case. The same phenomenon happens in 
the Einstein- Yang-Mills system where non-singular finite energy solutions have been discovered 
some time ago @ . 

The above YMH model without a dilaton field is known to possess nonsingular, finite energy 
solutions, describing magnetic monopoles Q. In the present work we study static, spherically 
symmetric solutions of the DYMH theory. Our results show a striking similarity to those ob- 
tained in the corresponding YMH theory coupled to gravity (EYMH) In the DYMH theory 
there is a finite energy abelian solution for all values of the dilaton coupling a. Based on nu- 
merical investigations there is a strong indication for the existence of a finite energy nonabelian 
monopole up to only a maximal value of a, am- The nonabelian monopole merges with the 
abelian one at a critical value of a, ac- In limit when the dilaton decouples the nonabelian 
solution joins smoothly the 't Hooft-Polyakov monopole. If the Higgs self coupling, /3, is in the 
interval [0,0.6] a^aiP) and adP) are different and with close analogy to the EYMH theory there 
are two different solutions if ac < a < am- 

In addition to the 'fundamental' monopole there is good numerical evidence that a countable 
family of globally regular solutions exits for < a < \/3/2 independently of /3. As their energy 
is higher they can be interpreted as excitations of the fundamental monopole. In the limit 
when the dilaton filed decouples their energy diverges. There is another limit when the DYMH 
theory reduces to the DYM model and then these excitations tend to the solutions discovered 
in refs. |,§. 

The action of our model is given by 

S = J d'^x !^^d^^d>^ip - e^'^^^F^.F-^^'' + ^{D^^r{D^^<^r - e-^^'^Vi^^^ , (1) 
where a = 1,2,3, and the Higgs potential is 

= -t;2)2_ (2) 

8 

The couplings of the scalar field ip has been chosen so that the action (||) possesses the following 
dilatational symmetry 

xf^^e^'^x'', A^^e-^'^A^, ip^^ + c, S ^ e^^^S , (3) 

which can be used to fix the value of ip arbitrarily at any point. With the 

$ u<I>, X —X, (f v(p, S -S, 

A — > 5^/3^, K -a 

V 

rescaling of variables the only remaining dimensionless parameters are 

^^Mw . Mh 
g Mw 

where Mw = vg and Mh = are the characteristic masses of the theory. Without further 
restrictions one can assume a > 0. Let us mention here that the action of the DYMH theory 
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(|l]) can be obtained from the Einstein- Yang-Mills-Higgs action by restricting the metric to the 
special conformstatic form ds"^ = e^'^dt^ — e~'^'^dx^dx^ . 

We use the minimal static spherically symmetric ansatz with zero 'electric' field 

Ag^O, A^ = eaij^iWir)-l) , (6a) 

^^ = H{r)—, (6b) 
r 

ip = ip{r). (6c) 
After the rescaling (Q) the action ([l|) reduces to 
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E = I dr{ irV' + e=^"^ + + W'H^ + ^ + i^r^e-^'^iH^ - if \ (7) 
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using the ansatz (|6|). Varying we obtain the equations of motion: 

^r'^')' = 2ae'^^{w" -^rh-'^^{H' - If , (8a) 

r^W" = W{W^ - 1 + r^e-^'^'PH^) - lar^^f'W , (8b) 
{r^H'Y = h(2W'^ + ^r2e-2"^(i/2 _ i)) . (8c) 

Solutions regular at the origin must satisfy the following boundary conditions (b.c): 

H = ar + 0{r'^), W = I - br"^ + 0{r^), 

i.e. there is a three parameter (a, 5, ipo) family of regular solutions at r = 0. The local existence 
can be proved following the procedure discussed in |^] 
For r — > cxD the corresponding 'regular' b.c. are 

ip = iPoo-^ + 0{^) , (10a) 

W = Be'^'P (l + 0{-)] , (10b) 



H = l--e~^'^p(l + 0(-)] for/3<2 (10c) 
r V r J 

H = l-^J^^^e~'^^{l + Oil)) for/?>2, (lOd) 

parametrized by {foo, d, B, C), where fj, = 6""*^°°, p = r + adlnr. For /? > 2 we cannot fully 
parametrize the stable manifold at r = oo. Exploiting the virial theorem and using Eqs. (|^), (^) 
and ( [lO| ) one obtains for the energy (|^) 

oo 

E=- [dr{r\')' = -{r\')^. (11) 
a J a 



We discuss first an exact two parameter family of solutions of Eqs. (^ which is going to 
play an important role in our further analysis, since it descibes the asymptotic behaviour of 
nonabelian solutions. Consider the singular abelian monopole, W = H = 1 then the general 
solution of (^) takes the form 

V.a = --ln|^sinhUi- + ^))| (12) 
3 
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Figure 1: The a dependence of the fundamental monopole solution for /3 = 



with A^,M G R. For M > and A real it is regular for all r > 0. Its energy, E Eq. (^), is 
infinite unless ^ = 0. For ^ = the total energy of the solution (|l^ ) is finite, E = Mja. When 
M < and A is real the solution is singular at r = —aM and the energy diverges. If A is 
imaginary for any value of M the solution becomes singular at some r > 0. 

Using the dilatational symmetry from now on we set ipoo = 0. Note that then the energy 
scales as ^ exp(— a^Joo)-^. The finite energy abelian solution for Lp^o = reads 

1 , / «N 
^= In 1 + - , 

a r 

and its energy is simply E = \la. It satisfies the following first order differential Eq. 

rV = ±e"^ , (13) 

which is derivable from a Bogomolny type bound ||5| in the background W = H = 1. 

The asymptotic behaviour of finite energy nonabelian solutions is described by (O). With 



the help of Eq. (11) we arrive at the following useful formula for the total energy of globally 



regular nonabelian solutions: 

E = - coshiA/M) , (14) 
a 

where A, M are parameters of (pa. describing their asymptotic behaviour. For globally regular 
solutions of {|8|) satisfying the b.c. @, (pX|) the energy (|^ is a monotonically increasing function 
of the Higgs coupling (3 and monotonically decreases with the dilaton coupling a. (We remark 
here that the condition ^poo = is essential to prove these facts.) 

We have numerically integrated a suitably desingularized version of Eqs. (^) from r = using 
the b.c. (^). We set ^po to zero in order to have to vary only two parameters. By adjusting the 
parameters a, h we suppressed the divergent modes in the asymptotic region according to (|l7|). 
The resulting solutions can then be transformed to satisfy (poo = by the dilatational symmetry 
®. 

We found that there exist a fundamental nonabelian monopole solution for all (3 below a 
maximal value of a, am which smoothly joins the 't Hooft-Polyakov monopole in the a — > 
limit. Some solution curves for /J = are shown in Fig. |l|. For the parameters see Table |l[ 

There is a critical a value, Oo where the dilaton field becomes logarithmically divergent 
while W and H tend to some nontrivial functions as can be seen on Fig. |l|. After shifting ip^o to 
the fundamental monopole tends to the finite energy abelian one as a ^ Oc. For < 0.06 an 
is seen to differ from (see Table ^) . It means that for ci given (X in the interval \cxc^ oixxi 

1 there 
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Table 1: Table 2: 



Parameters for the fundamental monopole solutions in Fig. 1 




(3-' 




am 


a 


a 


h 




a. ■ E 





1.39938 


1.4088 


0.0 


0.33333333 


0.16666666 










0.02 


1.37874 


1.3803 


0.1 


0.33324422 


0.16674419 


0.050094 


0.099853 




0.03 


1.36950 


1.3702 


0.2 


0.33297626 


0.16697893 


0.100758 


0.198827 




0.04 


1.36088 


1.3612 


0.5 


0.33107418 


0.16872021 


0.262673 


0.481732 




0.043 


1.35839 


1.3586 


0.7 


0.32884174 


0.17095758 


0.387902 


0.650049 




0.05 


1.35279 


1.3529 


0.8 


0.32741363 


0.17252229 


0.460048 


0.725582 




0.06 




1.3452 


1.0 


0.32388447 


0.17699856 


0.638684 


0.855180 




0.07 


"m 


1.338 


1.1 


0.32176986 


0.18029088 


0.758343 


0.907479 




0.08 




1.33 


1.2 


0.31944319 


0.18485325 


0.919110 


0.950109 




0.1 




1.32 


1.3 


0.31705890 


0.19195967 


1.172960 


0.981707 




0.15 




1.29 


1.4 


0.31657705 


0.20987214 


1.950448 


0.999523 




0.2 




1.27 


1.4088 


0.31898696 


0.21924215 


2.582925 


1.000092 




1 




1.11 


1.4 


0.32380520 


0.22987901 


5.289654 


1.000000 
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0.98 


1.39939 


0.32402298 


0.23025275 


8.461582 


1.000000 




oo 




0.89 



are two different nonabelian solutions with different energies. The function am(/3) decreases 
with increasing /? from ani(O) ~ 1.4088 to am(oo) ~ 0.89. 

There also seems to exist a countable family of globally regular monopole solutions indexed 
by the number zeros of W{r) for all /? and < a < V3/2. They can be interpreted as radial 
excitations of the fundamental monopole. In the a ^ limit after a suitable rescaling they 
can be identified with the previously discovered DYM solutions @],|^]. We illustrate some of 
these excited solutions with one and two zeros for a various values on Figs. § and |3|, while the 
corresponding parameters are listed in Table ^ and Table 



Table 3: 



Parameters for the first excited monopole solution in Fig. |2 




a 


a 






a ■ E 


DYM 




0.26083015 


1.711412 


0.8038078 


0.05 


0.66718265 


0.26286356 


1.731264 


0.8507008 


0.1 


0.77749002 


0.27004539 


1.803886 


0.8904364 


0.2 


0.99147112 


0.30099818 


2.170909 


0.9448188 


0.3 


1.00351436 


0.33435729 


2.730953 


0.9720677 


0.5 


0.77902186 


0.37296818 


3.994232 


0.9932192 


0.7 


0.59081576 


0.39281185 


5.919563 


0.9992733 


0.8 


0.51698205 


0.39916099 


8.207917 


0.9999547 


0.866 


0.47164657 


0.40215328 


26.31413 


1.0000000 



When a — > \/3/2 for the excited solutions the dilaton diverges logarithmically again while 
the zeros of W go rapidly to infinity, except for the innermost one. For \/3/2 < a < 1 there is 
a surviving solution with a single zero of W and with divergent ip existing up to a = 1, where 
W develops an extra divergent mode, so for a > 1 this solution is not expected to exist. There 
also exists another type of limiting solution when the number of zeros of W goes to infinity for 
all < a < \/3/2 and the dilaton field diverges logarithmically. 

If one shifts, however (/?oo to then all of the excited monopoles merge with the finite energy 
abelian solution for a — > \/3/2- The fundamental and the excited monopoles up to six zeros are 
plotted on Fig. ^. In order to better display the higher zeros of W we plotted ^/rW . Notice 
that the newer zeros appear nearer and nearer to the origin. 
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Figure 2: The a dependence of the first excited monopole solution for /3 = 




10 

ln(l + r/a) 



10 



DYM 




0.05 




0.1 




0.2 - 




0.5 




0.840 

























10 12 

ln(l + r/a) 



Figure 3: The a dependence of the second excited monopole solution for /3 = 
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Figure 4: The fundamental and the first six excited monopole solutions for a = 0.2 and /? = 
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Table 4: 



Parameters for the second excited monopole 


solution in Fig. ^ 


a 


a 






a ■ E 


DYM 


— 


0.35351801 


3.373903 


0.96559852 


0.05 


6.0251373 


0.36759603 


4.169777 


0.99251583 


0.1 


3.9102367 


0.37450939 


4.996143 


0.99673418 


0.2 


2.0946637 


0.37885368 


5.999793 


0.99882187 


0.3 


1.4148625 


0.38161925 


6.815222 


0.99949621 


0.5 


0.85071037 


0.38762800 


8.666195 


0.99993064 


0.7 


0.60036405 


0.39522735 


12.25580 


0.99999862 


0.8 


0.51818563 


0.39944391 


17.77525 


1.00000000 


0.84 


0.48948129 


0.40111050 


25.77133 


1.00000000 



The structure of the solutions can be understood from the linearization of the field equations 
around the abelian monopole (|l^). Using the 

W = w{r), H = l + ^, if = ip,{r) + ip , (15) 

r 

variables, from the linearized field equations one finds for l/M <^ a/r <^ 1/|^| (where ip^ — 
ln(r/a)/a) that the solutions are well approximated by: tp ~ e'^'^'^ , h ~ e^'^'^, w ~ e"^™"^ 
wherer = Inr and the 'frequencies' are 

A^ = (l;-2), Xh = ^{l±^Jl + Aa^(3^) , \^ = ^Ja^ -3/A. (16) 

So we see that in this 'middle' region W oscillates with an amplitude decaying like l/\/r. If 
M — > oo (i.e. aE 1) this region streches out to infinity, while W has more and more zeros when 
a < \/3/2. This behaviour is a similar to the one found in the DYM case in the interval where if 
grows logarithmically. If a > 1 1^ we do not expect the corresponding solution to exist. In the 
asymptotic region defined by r ^ |a74coth(A/M)|, r ^ (where /U = | sinh(j4/M)/^| and 
~ — i In/i) the linear corrections /? 7^ are characterized by ~ e^^'^ , w ~ e^"'^, h ~ e^'^'^, 
with 

A^ = (0;-1), Xh = ±fi\P\, A^ = ±/i, (17) 

while h ~ e^ft^ and Xh = (0; +1) for /3 = 0. 

The energy of the solutions goes rapidly to 1 if the number of oscillations of W increases 
(see Table |5|). We have determined the energy of the solution by fitting the parameters A, 



Table 5: 



Parameters of the solutions for a = 0.2 and /? = in Fig. |^ 


n 


a 


b 




a ■ E 





0.33297626 


0.16697893 


0.1007578 


0.1988268 


1 


0.99147112 


7.52495447 


10.854544 


0.9448188 


2 


2.09466373 


9.47134195 


29.998964 


0.9988219 


3 


2.14396020 


9.51995711 


48.727339 


0.9999720 


4 


2.14545164 


9.52111922 


67.372640 


0.9999993 


5 


2.14549487 


9.52114737 


86.014682 


1.0000000 


6 


2.14549609 


9.52114805 


104.6566 


1.0000000 



M in Eq. (12) in the asymptotic region using formula ([l^). The energy determined this way 
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Figure 5: a dependence of the energy of the fundamental, the first, and the second excited 
monopole solutions, and its detailed view for the fundamental monopole near for /3 = 
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Figure 6: Detailed view of the a dependence of the energy of the fundamental monopole 
solution near for = 0.02 and for = 0.043 



contains only exponentially small corrections. We have also plotted the a dependence of the 
energy, i?(a), (rescaled to ^poo = 0) on Fig. ^ for /3 = 0. The similarity of Fig. ^ to Fig 3 in 
Ref. Q where the masses of gravitating monopoles are plotted as a function of the 'gravitational 
coupling strength', Mw/Afpianck, is indeed striking. On Fig. ^£'(0;) for = 0.02 and 0.043 is 
shown. 

For not too large /3 values (/? < 3) aE of the fundamental monopole becomes larger than 1 
unlike for the excited ones. For (3 = 00 this maximum is 1. 

We make finally some remarks on the stability of the solutions. The 't Hooft-Polyakov 
monopoles are stable since they are solutions with minimal energy |1C]. It is natural to expect 
the fundamental monopoles to remain stable in the DYMH case for a > 0. For sufficiently small 
P, however, where the mass of the fundamental monopole is larger than that of the abelian one, 
the nonabelian solution is expected to become unstable against large perturbations. If Uc / Om 
there is a bifurcation point where the linear stability of the solutions can change. In the EYMH 
case this change of stability has been shown in The excited monopoles are expected to be 
unstable for all a since their energy is significantly larger than that of the fundamental ones. 
This heuristic argument is strengthened by the fact that in the a ^ limit their counterparts 
are known to be unstable |5|, ^. 
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